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We propose a model for two (d + l)-dimensional directed polymers subjected to 
a mutual ^-function interaction with a random coupling constant, and present an 
exact renormalization group study for this system. The exact /3-function, evaluated 
| through an e(= 1 — d) expansion for second and third moments of the partition 



function, exhibits the marginal relevance of the disorder at d = 1, and the presence 
of a phase transition from a weak to strong disorder regime for d > 1. The lengthscale 
exponent for the critical point is v = 1/2 | e |. We give details of the renormalization. 



We show that higher moments do not require any new interaction, and hence the /? 



function remains the same for all moments. The method is extended to multicritical 
systems involving an m chain interaction. The corresponding disorder induced phase 



transition for d > d m = l/(m — 1) has the critical exponent v m = [2d(m — 1) — 2] . 
For both the cases, an essential singularity appears for the lengthscale right at the 
upper critical dimension d m . We also discuss the strange behavior of an annealed 
system with more than two chains with pairwise random interactions among each 
other. 
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I. INTRODUCTION 



Attempts to study the effects of randomness, especially if one requires averages of ther- 
modynamic quantities, have led to many new techniques, concepts, and, probably, contro- 
versies. In order to get a clear idea about random systems, in recent years, a directed 
polymer (DP) in a random medium seems to have emerged as the consensus candidate for 
the "simplest" random model. [0,fM,l5]|5[] We here propose a still simpler problem of DPs with 
random interaction that can be solved using an exact field theoretic renormalization group 
(RG) approach. |7J This, we believe, is highly significant since RG is the general framework 
to study and to understand, through the fixed point spectrum, the universal aspects of any 
model. 

Directed polymers in (d + l)-dimensions are random walks directed along a particular 
direction, say z, with fluctuations in the transverse <i-dimensional space. DPs are of con- 
siderable interest and have attracted a lot of attention as a simple statistical mechanical 
model because of its relevance and applicability in unifying a wide variety of seemingly dis- 
parate systems. These include the flux lattice melting problem in high T c superconductors, 
H commensurate-incommensurate transitions, wetting transition, [|K],[ll|] vertex models, 
fT2|j polymeric nematics, biomembrane phase transitions, |L4| interface growth, [H| etc. 
Many problems of conventional polymers [self avoiding walks (SAW)] like collapse, adsorp- 



tion etc, have exactly solvable counterpart in DPs. |16| The RG analysis of a pure system 
of interacting DPs gives enough insights through the evaluation of the exact /^-function to 
all orders in perturbation series. HT7j|T2"l , |l8f These systems of DPs with pure short range in- 
teractions are almost completely solved, and, for example, have led to several exact results 
for vertex models. |12| . |T9| 

There are many efforts and activities in the field of polymers with random interaction 
20| or in random media [21]. The analogous DP problems are expected to be simpler. 



For example, a DP in a random medium, which through a nonlinear mapping describes 
many aspects of interface growth, has been studied upto one loop in the momentum shell 
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technique. |15| There were also attempts for solving the many chain system in a random 
medium in the context of high T c superconductivity. |22] It is the directedness that helps 
in setting up the DP problem, both analytically and numerically, as oppposed to the SAW 
problem in random media. Several results for the DP problem are known in general, though 
exact or rigorous results are rather few. || Apart from these random media problems, the 
other category of problems involve polymers with random interactions in the context of, say, 
disordered heteropolymers. [^(J Here again, a DP with random interaction turns out to be 
simpler. |23[ 

Our model @ has similarity with the second category of problems. It deals with a random 
mutual interaction among the chains with the randomness in the coupling constant of the 
interaction. The randomness is only along the length of the chains and does not depend 
on the transverse d dimensional coordinates. The specific charateristic of the randomness 
as well as the directed nature of the polymers enable us to solve the model exactly. We, 
furthermore, show that this model, inspite of its simplicity, captures many of the essential 
features such as marginal relevance, existence of a disorder induced phase transition, etc. as 
known, e.g., for the interface growth problem, DP and SAW in random environments etc. 
PPT[] The two dimensional wetting phenomenon is also analogous to our proposed system 



il| . [2^| - though our model (and the solution) is for general d. 

We define the model in the next section, and to put things in the proper context, the 
aim and the outline of the paper are given there. 



II. MODEL 

From the definition of DPs, it follows that a projection of a DP in the transverse d- 
dimensional space is an ordinary polymer with z representing the contour variable which 
is equivalent to the steplength in a discrete case. In the path integral formulation the 
dimensionless hamiltonian for two such DPs, each of length N, interacting through a random 
mutual shortrange interaction can be written as [|7| 
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where r^z) is the (^-dimensional position vector of a point of chain i at a contour length z, 
and ri2(-z) = ri(-z) — ^(.s). The first term that comes from the chain connectivity is the 
entropic contribution, and corresponds to free chains. |25] The second term is the two chain 
interaction at the same chainlength through a short range potential V^(r). We introduce 
the randomness through the coupling constant. It has a pure part Vq and a random part 
Vob(z) which varies only with z (the length along the chain). It is chosen in this way so that 
b(z) is dimensionless. At this stage, for generality, we keep V as a short range potential. 
Later on, specific calculations would be done with a 5-function potential. Also, starting 
with a short range potential has certain mathematical advantages like avoiding powers of 
distributions, as we will see below. One can also think of this problem as a nonrelativistic 
quantum problem of particles with time (z) dependent interaction potential - a description 
we do not find very illuminating. 

One of the simplest but nontrivial choices for the distribution of the randomness is a 
Gaussian one: 

P(b(z)) = (2ttA)- 1/2 exp[-&(z) 2 /(2A)], (2.2a) 
(b(z)) = 0, and (b(z 1 )b(z 2 )) = A 5( Zl - z 2 ). (2.2b) 

Here the randomness is uncorrelated in nature and is described by the variance A. Choosing 
a zero mean for b(z) is not a restriction because any nonzero (b(z)) could be "gauged away" 
by absorbing it in the pure part. 

So far we have discussed only two body interactions. For DPs it is known that even 
pure many body interactions, representing special multicritical points, can also be handled 
exactly. [f^,|19|,[n| It is possible to study the disordered versions of these multicritical systems. 
The hamiltonian for the m-th order multicritical point, involving only m-body 5-function 
interaction, is 

1 r N m fdr (z)\ 2 r N m ~ 1 
H ™=2J dz ]^\^) + ^ dzv ™{ 1 + b ^ II<K r -+i(*)) (2-3) 



2 N , 2- 

dri(z) \ ( dr 2 (z] 



r-N 



where, as before, b(z) is the random part. We come back to this multicritical situation in 
section [VT]. It is also possible to define more general systems by putting the lower order in- 
teractions in Eq. |2.3| with independent random coupling constants. Such a hamiltonian can, 
in principle, describe the approach to the multicritical points. However, such complicated 
cases are not discussed here. 

One possibility of getting a random interaction for, say, the two chain case, is to take 
"charged" DPs, with random charges qi{z) for the ith chain, vob(z) = qi(z)q2(z), and in- 
teractions of charges only at same z. If the charges are in thermal equilibrium with the 
polymers, a simple quadratic hamiltonian for the charges can be taken as proportional to 
/ b(z) 2 dz, and it is to be added to the hamiltonian of Eq. [2.1| . The partition function one 
gets from this full hamiltonian is really equivalent to (Z) (annealed case) evaluated with the 
hamiltonian of Eq. [2.1| but averaged over the distribution of Eq. [2.24 The more complex 
situation is the quenched average which requires, e. g., the average of the free energy. 

One of the standard approaches for random systems is to proceed through the evaluation 
of the quenched free energy using the replica trick [|TJ 

[Z n ) - 1 



(InZ) = lim 



n^oo 



n 



which requires evaluation of (Z n ). The importance of these moments can be realized through 
the following expansion 

oo 

<lnZ> = ln(Z> + £ n-\[Z/(Z)-l\ n ). 

n=2 

Such an expansion makes sense if and only if the various cumulants of the partition function, 
with respect to the disorder distribution, do not grow too rapidly with n. In such a case 
there will not be much qualitative difference between the quenched and annealed cases. 
This, in turn, suggests that to look into the possible differences one can study the various 
moments of the partition function. In addition, the moments can be looked upon as the 
charateristic function for the probability distribution of In Z. [[^J They are, therefore, of 



interest in themselves. [£8] This is the approach we take in this paper. Our analysis is not 



yet enough for the analytic continuation in n to < n < 1 regime, as one would need for the 
free energy. This is not a deterrence as important information can be gathered even from 
the integral moments. 

Instead of evaluating the quenched free energy, we systematically study the behavior of 
(Z), the second and the third cumulants. The first moment describes the behavior of an 
annealed system while the higher cumulants would show the nature of fluctuations. For each 
case, we do the averaging exactly before the configuration sum (or "path integrals") to define 
an effective hamiltonian for that particular cumulant. This "pure" effective hamiltonian 
is then treated by perturbative renormalization. The effect of disorder is felt through the 
generation of new terms in the effective hamiltonian. The RG analysis helps in examining the 
flow of these terms as the lengthscale is increased, thereby, showing the marginal relevance 
or irrelevance of the disorder. 

In section [jT| (Z) is discussed, while (Z 2 ) and (Z 3 ) are done in section |TV| and [V]. 
Though the derivation of effective hamiltonians preceeds perturbative analysis in sections 
[IT]] , [TV] and [V], it is instructive to start with the original hamiltonian, do the perturbation 
analysis and then do a term by term disorder averaging. Such a procedure not only shows 
how the new terms are generated but also acts as a cross check. This is discussed in Appendix 



A. Appendix B discusses many of the details needed in section [TV]. The random multicritical 



case is discussed in section VI. The annealed case of three and four chains is discussed in 



section VII. The paper ends with a discussion and summary in section VIII. 



III. (Z) 

We show in this section that the annealed case can be reduced to a pure problem. We 
add that this reduction is special for two chains. Had we started with more than two chains, 



say three or four, with the same random pairwise interaction as in Eq. |2.1| , the annealed 
case would be completely different from the corresponding pure case, and would have a much 
richer structure. A particular case is discussed in section |V11. 
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The partition function, in the continuum approach, for a system of two chains, given by 
the hamiltonian of Eq. |2.1|, is 



Z = J Dri Dr 2 exp(— H) 

where / T)r\Dr% stands for the sum over all configurations of the two chains. A straight- 
forward averaging of Z using the probability distribution of Eq. |2.2a| defines an effective 
hamiltonian 7i e s such that 

(Z) = J D n Dr 2 exp(-Kff), (3.1) 

and it is given by 

2 

i r N 2 fdr(z)\ r N v 2 A r N 

n ^ = U ^M^r) +Vo l dzV{ri2{z)) ~ 2 L dzV2 ^ z ))- ( 3 - 2 ) 

It appears from the above expression of the effective hamiltonian that an attraction is 
generated between the two chains. We find it instructive to follow another approach of 
perturbation expansion of the interaction term starting with the original hamiltonian |2.1| . 
This helps us in visualizing the origin of the disorder induced attraction. This is done in 
Appendix A. 

Now, since any short range potential under renormalization maps onto a 5 function 
potential, we can take the "minimal" effective hamiltonian for (Z) as 



rN 

2 Jo 



jy z !^je^ ^ j:dz5Mz)Y (3 . 3) 



where Vq is the reduced coupling constant which takes care of the attraction described earlier. 
We beleive that the large length scale properties as described by Eq. [T3] is same as that of 
Eq. p.2| . If necessary, we can restrict the strength of the disorder so that Vq, which represents 
the effective coupling between the two chains, is positive (i.e. repulsive interaction). Now 
the problem reduces to a relatively simple situation where the two chains interact with a 
pure 5-function interaction with a reduced coupling constant v . The solution of this pure 
problem is known, and is used below. |T7|Jil 



For the sake of completeness we quote the relevant results from from Ref. [fTTf . The 
perturbation series for the connected part of the annealed partition function (Z) c to all 
orders in vq is 



(Z) c = NVv 



£i (47r)" d / 2 r(2 + ne') 



(3.4) 



where V is the d- dimensional transverse volume, and r(«) is the standard gamma function. 
The exact /3-function for the renormalized coupling constant u (with uq = voL 2 ~ d as the 
bare dimensionless coupling constant) 

(3{u) = L% = 2e'u ( 1 - -^-) . (3.5) 



dL V 4vre' 

Note, here 2e' = (2 — d) replaces e of Ref to avoid later conflict of notation. 



The flow diagram for the dimensionless coupling constant u is shown in Fig. |l|. The fact 
that for d < 2 any small attractive interaction is able to form a bound state is reflected by 
the flow to the nonperturbative regime for any negative u. The repulsive or the positive u 
region is dominated by the stable fixed point u*{— Ante'). For d > 2 there exists a nontrivial 
unstable fixed point u = u* which seperates the bound and the unbound states for the two 
polymers. In short, the unstable fixed point represents the critical point for the binding- 



unbinding transition. The exponents are known and can be found in Ref. |32| . For example 



for 1 < d < 4, the length scale exponent is 1/| e' \. 



IV. {Z 2 ) 

The evaluation of (Z 2 ) closely parallels that of the previous section. However, unlike the 
(Z) case, new terms are generated here in the effective hamiltonian. An RG analysis is done 
to get a detailed account of the effects of these new terms. 

A. Effective hamiltonian 



The averaging for (Z 2 ) with the hamiltonian in Eq. |2TI| needs completion of the per- 
fect square associated with b(z). As in the replica analysis |l| where one needs n replicas 
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("copies") of the original system in evaluating (Z n ), we require four chains for {Z 2 ), a 
pair {3,4} as a replica of the original pair of chains {1,2}. Therefore, we write, restricting 
ourselves to 5-function potentials, 

(Z 2 ) = J {[Dn exp(-W 2>2 ), (4.1a) 
i=i 

where 

H 2 , 2 = H + H 1 + H 2 , (4.1b) 

with 

*-5jr*g(^)" 



denoting the four chain free part, and 

i-N 



H l= v dz [5(r l2 (z)) + 5(r u (z))} (4.1d) 
Jo 

and 

H 2 = -f / dz 6(r 12 (z)) 5(r u (z)) (4.1e) 
Jo 

representing the interactions among the chains, with fo = v^A. 

In Eq. |4.1b| , Hi denotes the repulsive interaction between the chains of a particular 
pair ("intra replica", {12} and {34}, no cross coupling) at the same chain length, and is 
identical to the interaction term used for (Z) as discussed in section [I II] and Appendix A. 
The other term H 2 couples the two pairs of chains {12} and {34} ("inter replica" term), 
and is the crucial term for our analysis. Eventhough this is a four chain interaction, it is 
distinct from the multicritical type interaction of Eq. |2.3j . It cannot be interpreted directly 
as a conventional interaction between the two pairs. Rather, there is a lowering of "energy" 
of the system if the partners of each pair {12} and {34} meet simultaneously at the same 
chain length but not necessarily at the same point in space. This can also be interpreted as 
a special correlation so that an encounter of {12} at a chain length z favors an encounter 
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for {34} right at the same length z. A tendency to achieve this kind of configurations leads 
to all the nontrivial effects of the disorder. 

The coupling constant of H 2 , fo in Eq. [4.14 appears to be similar to that of the attractive 
interaction which is present in Hi, Eqs. [4. ldj and |3.2| , but they require separate treatments 
in the RG analysis (see below). As discussed in Appendix A, the term proportional to VqA 
in vq is reduced by a cutoff volume factor Q needed to define the 5 2 term properly. Because 
of this reduction VqA/Q differs from fo in dimensionality and matches properly with vq, the 
coupling constant of the starting S function interaction of the two chains. 

The standard dimensional analysis for dimensionless hamiltonian shows that [vq] = 
L d ~ 2 , [A] = L 2 and hence [v^A] = L 2d ~ 2 where L has the dimension of length. There- 
fore, the upper critical dimension for fo is d = 1 which also appears as a special dimension 
through the divergences in the e(= 1 — d) -expansion to be discussed below. From this simple 
dimensional analysis it also follows that the coupling in H 2 differs from that of Hi , as already 
mentioned. The special dimensionality d — 1 which is associated with H 2 is more important 
in the context of fluctuations in the partition function. In the absence of this term there is 
no special effect of disorder, which, in turn, also means that the quenched and annealed free 
energies would become equal. 



B. PERTURBATION SERIES 

To study the effect of H 2 , we develop a perturbation series for (Z 2 ) in f . The divergences 
that appear are absorbed by renormalization through an e-expansion. We show that this 
renormalization can be carried out exactly to all orders. For simplicity, this is done first for 
the Vq — case. We then show that these divergences at d = 1, arising only due to H 2 , 
remain unchanged even if we include Hi i.e when fo 7^ 0. The (3 function evaluated exactly 
to all orders in perturbation series and other essential features are identical for both f = 
and fo 7^ cases. 
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1. v = 



Let us consider first, for simplicity, the case when = 0. This means that there is no 
mutual two chain interaction. We consider only the connected part {Z 2 ) c = (Z 2 ) — (Z) 2 , 
the second cumulant of the partition function. As in the previous section, the calculation 
can be done in the real (chain) space. But, at this point, we prefer the Laplace space 
(Laplace transform with respect to the chain length) because it is advantageous for later 
considerations especially with vq ^ 0. We define 

Z= / dNe- sN (Z 2 ) c (4.2) 
Jo 

the Laplace transform of {Z 2 ) c with respect to the chainlength N, the Laplace conjugate 
variable being s. 

The loops in the perturbation expansion are shown in Fig ||]a upto third order in the 
interaction. The individual pairs of chains are represented by thick lines. The horizontal 
wiggly lines in these diagrams stand for f . Such a representation is possible because the S 
function in H2, Eq. [4.14 forces the members of a pair to have the same r, z coordinates. 



Each chain is described by the free distribution ("propagator") G(vf — | Zf — Zi) = 
[2ir(zf — Zi)]~ d l 2 exp[— (rj — v i ) 2 /2( y Zf — Zi)\ with end points (rj, zf) and (r^, zf). Two chains 
are therefore described by 

G 2 (r I z) = (Aixz)- d/2 G(r | z/2). (4.3) 

This G 2 is the propagator for the thick lines. At each wiggly line, connecting four chains (all 
four having the same chainlength z), there are two integrations over the spatial coordinates 
of the two separate pairs of chains (thick lines). The loops formed out of the wiggly lines 
are only responsible for the divergences at d — 1. 

In order to trace the algebraic origin of the singularity, note that, by very nature of the 
interaction, the spatial integrations associated with the two thick lines are independent of 
each other. Each section of the thick lines, with Zi, z 2 as the end points, in a loop formed with 
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the wiggly lines, contributes (z\ — z 2 )~ d / 2 from the identity in Eq. |4.3| . Since the interaction 
demands same z for the two thick lines, the z integrals involve (z\ — z 2 )~ d type factors whose 
Laplace transform would contribute T(l—d) with pole at d — 1. The two independent spatial 
coordinates which are left out after the successive use of the normalization / dr G(r | z) — 1, 
lead to a V 2 factor for each diagram. The convolution nature of the z integrals, thanks to 
the time ordering, leads to a simple product of the individual Laplace transforms of the 
integrands, resulting in a geometric series for Z. The details of the evaluation of a few 
diagrams are given in Appendix A. 

Defining the dimensionless coupling constant ro through an arbitrary length scale L as 
r = f I/ 2e (47r)~ d , e = 1 — d, we write the series for Z to all orders in r as 

oo 

Z |, 0=0 = (4vr) d V 2 s- 2 [r + £ r " +1 {sL 2 )^ T n (e)]. (4.4) 

n=l 

It is clear from the above expression that there is a divergence at d — 1 at each order (> 1). 
This is tackled by renormalization below. 

2. v + 

In the above analysis, we have taken vq = 0. We now include vq and show that the 
singularity structure around d — 1, as in Eq. [4.4j , remains unaffected. 

In this case there are both intra-pair and inter-pair interactions. The intra pair interac- 
tions, i.e., the mutual short range ^-function interactions among the members of the pairs 
are represented by dots on the thick lines. The basic idea of the procedure adopted is to 
show that the dots can be absorbed by dressing the "propagators". The original propaga- 
tors (thick lines of = case) G 2 is modified by the dots but not trivially. The "dressing" 
factor depends on whether the chains are open or tied at the ends. One therefore needs two 
types of dressed propagators, Gm for the thick lines in the loops, and Go for the same in 
the outer legs not involved in the loops. (See Fig c and d.) Since these involve only two 
chain interactions, the singularities are at d = 2 of the type r(l — d/2) as in Eq. |3^4] . These 
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dressed propagators are to be used, as appropriate, in the skeleton diagrams of the vq = 
case without the dots (see Fig |2] a). We just qoute the forms of these propagators below - 
the details can be found in Appendix B. 

For the one in which the two participating members of a thick line are tied together at 
both the ends [(o, 0)and(r, z)] of the line, the form of the dressed propagator with n meetings 
(dots) is [Fig. |d ] 

Gt\r | z) = (-r;o) w (47r)-^ 1 ) rf / 2 r( i (w+ ^ ) J e/) g ^- 1 G(r | z/2). (4.5) 
There is translational invariance in both r and z. With e' = (2 — d)/2, this form surely 



reduces to Eq. [4.3| for n = 0. Similarly, for the other type in which the two members of a 
thick line are tied together only at one of the ends, integrations for the open end coordinates 
need to be done. The resulting dressed propagator for chains of length z with n intermediate 
dots (Fig 0c) has the following form 

G (n) (z) h^Wg U 6) 

G ° [Z) ~ r(l + ne')(47r)^- 

With n = 0, Gq\z) = 1, as it should be, by the normalization of the distribution function 
G(t I z). Also G { o\z) has no space dependence. 

In a diagram of a particular order in fo, the thick lines can have arbitrary order n in vq 
(i.e. arbitrary number of dots). All such diagrams differing only in orders of Vq are combined 
together by summing over n. The full dressed propagators are 

G M (r | z) = E^m ( r I *), and G Q (z) =J2G%\z). (4.7) 

n n 

Unlike the ^0 = situation, these two propagators replace the inner and the outer thick 
lines respectively. The subsequent procedure is almost similar to the previous case, including 
the origin of the V 2 factor, and the use of the convolution theorem in the Laplace space. 
The series for Z is given by 



Z |^o= V 2 G o {s)r 



1 + E (r G M (s)) n 



n=l 



Qo{s) (4.8) 
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where Qo(s) = CGq(z), Qm{s) = CG 2 M (z), (£ being the Laplace transform with resprect 
to z) and Gm{z) = J drGM^ \ z). The "same z" requirement of the 5 function of H 2 
combines the propagators of the thick lines. Hence G 2 Q and G 2 M . The two outer pairs of legs 
contribute the two Qo{ s ) factors. In terms of the dimensionless coupling constant r , and 
Uo(— vo(4:Tc)~ d / 2 T(2e')s~ e '), the above expression becomes 



Z |^ = (4tt) d V 2 s- 2 L- 2 *S (s) 



ro + J2ro +1 (sL 2 )-^(S M (s)r 



71=1 



So(s) 



(4.9) 



where 



Ms) - E r n K '" + 11 
w 1 1 



+ n 1 e')r(l + n 2 e / ) 



(4.10) 



and 



? , x _ r[(rai + w 2 )e / + 1 - d] 
M[S) ~^ } T((ni + l)e'n(n 2 + l)e') 



-u ) ni+n2 . 



(4.11) 



Details can be found in Appendix B. The reason it is written in the above form is that Sq 
and Sm start with 1 for n\ = n 2 = 0, to agree with Eq. P . It follows that the leading 
divergences at d = 1 in each order of Eq. [4.9| come from the n\ = n 2 = ... = term of 
Eq. Km and Eq. KTL 



C. Renormalization and New Criticality 



The divergences at e = in the series of Eq. O can be absorbed by the standard 
renormalization procedure. [30] In general, a renormalization through minimal subtraction 
would require absorption of the poles in e through 



r = r(l + a x r + a 2 r 2 + ...). 



(4.12) 



with a n = J2p=i 0"n,p£~ p and r as the renormalized coupling constant. In such a scheme, 
a n,p(p 7^ w ) terms are required to take care of the sub leading divergences. 

The formal similarity of the leading pole structure of Eqs [4.4| and [O] with that of 



Eq. |3]4] enables us to follow Ref. f[7] , [l2] yielding a p = (— e) p . The geometric series of Eqs. 



14 



]9| and |4.4| guarantees that the removal of the leading poles is sufficient to remove the 
subleading ones. The presence of the dots (Fig |2|b) through ra 1 ,n 2 7^ in Eq. [O] is felt 
through the changes in the subleading divergences. This does not pose a problem and can, 



indeed, be checked explicitly. Note that Sm of Eq. |4.11| has an expansion of the form 



S M = - + A + J2A p e p . (4.13) 
e P =i 

Taking a p = (— e)~ p , as needed to remove the leading poles, one can verify explicitly that all 
the poles are removed order by order, and the result does not depend on the explicit values 
of Aq, A\ etc. 

The (3 function is therefore exact to all orders in perturbation series and is given by 

P(r) = L— = 2(er + r 2 ). (4.14) 

There are two fixed points: (i) r = and (ii) r* = — e. The flows are shown in Fig |3|. 
The bare coupling constant r which originates from VqA, where A, the variance of the 
distribution, is strictly positive, requires a positive r. Therefore, the nontrivial fixed point 
for d < 1 in negative r is unphysical. It however moves to the physical domain for d > 1. 
See Fig. 3c. 

Exactly at d — 1, e = 0, r grows with length L as 



r(L) = r(0) 



l + 2r(0)ln^ 

Lj 



(4.15) 



r(0) being the coupling at length L . Hence, the disorder is marginally relevant, in agreement 
with Ref [24]]. For d > 1, there exists an unstable nontrivial fixed point at r =| e | which 



separates two distinct regimes of disorder. If we start with a strong enough disorder, on the 
right side of the fixed point, it increases with length scale, going beyond the perturbative 
regime. This is the strong disorder phase. On the other hand, the left side of the fixed point 
is the weak disorder regime, since r flows to zero (the stable fixed point). The unstable 
fixed point, therefore, represents a critical point - a novel phase transition induced by the 
disorder. 
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One way of achieving the abovementioned critical behavior is to change the strength of 
the disorder by controlling the temperature. The "strong disorder" phase ((InZ) 7^ ln(Z)) 
would correspond to the low temperature phase while the "weak disorder" phase ((InZ) ~ 
\n(Z)) is the high temperature one. The details of the critical behavior can be obtained by 
integrating the (3 function, 

2 I e 



1 



r(0)- 



r(0) 



L 



(4.16) 



For a small starting deviation AT = T — T c = r(0)— | e |, there is a lengthscale L ~ 
(- AT )-i/2|e| at which 

r in Eq. [4.16| diverges. This we can identify as a lengthscale £ 
associated with the critical point with the lengthscale exponent 



"=(2 I e |) 



-1 



(4.17) 



The divergence at e = is consistent with the essential singularity that follows from Eq. |4.15 



£~exp[l/(2AT)]. 



(4.18) 



A complete description of the critical point would involve an evaluation of various macro- 
scopic or thermodynamic properties. These would require a replica type analysis. It is 
tempting to believe that the correlation induced by H 2 in Eq. |4.1ej in the replica space 
distinguishes the two phases. We wish to come back to such replica analysis elsewhere. 



V. (Z*) 

The evaluation of (Z 3 ) leads to a six chain problem where, as before, an interaction 
involving four chains is generated which is attractive in nature. The effective hamiltonian, 
apart from the free part for six chains and mutual 5 function interaction, contains the 
following attractive terms (see Eq. [4.1 bp 

rN 

-f„ / dz [5(r 12 (z)) 5(r u (z)) + 5{r u (z)) S(r 56 {z)) + 5(r 12 (z)) a(r 56 (*))]. 
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Instead of (Z 3 ), we analyze the third cumulant (Z 3 ) c involving only six chain connected 
diagrams. Fig [| shows such diagrams upto fourth order in f . Contributions of these 
diagrams can be found out following the rules discussed in the context of (Z 2 ) in section 
P"V| . As an example, we give an explicit evaluation of Fig |]c which is 

f fN rzi rZ2 rz^ 

f / / dzi / dz 2 / dz 3 / dzi G 2 {ri 2 \ Z\ 2 ) G 2 {j\ 2 I z 12 ) X 

J{r,r',r"} JO JO JO JO 

G 2 { r> 23 I z 23)G 2 (r' 34: I 2 34 ) G 2 (r" 34 | z 34 ) 

In the above equation r, r', r" with appropriate subscripts denote the set of (^-dimensional 
coordinates for the three thick lines and J7 r r , r „| corresponds to the integrations over all 
spatial coordinates. As before, each thick line between two end points (r^, Zi) and (rj,Zj) is 
represented by G 2 {rij \ Zij) with Zij = Zi — Zj. The spatial integrations simplifies the above 



expression to 

r 4 (47r)- M / 2 V 3 f N dz 1 [ Zl dz 2 [ Z2 dz 3 P dz 4 z u d z 2 *' 2 z^ 

JO JO JO JO 



The series for the Laplace transform of {Z 3 ) c is given by 

[C(Z 3 ) C ] U 0=0 = V\Aitf d l 2 s- l - 3d l 2 T(e')[2r, 2 (sL 2 y 2 * + 4r 3 ( S L 2 )- 3 T(e) + 

6r 4 ( S L 2 )- 4 T 2 (e) + ...] (5.2) 
where, as before, s is the Laplace conjugate to the chainlength iV and r is the dimensionless 



coupling constant as defined before Eq. [4.4| . This series requires the standard renormalization 



procedure for removal of divergence at d — 1. Defining the renormalized r via 

r = r(l + air + a 2 r 2 + ...) (5.3a) 

it is found that 

a p = {- l/ef (5.3b) 
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absorbs the divergence at d = 1. 

It is interesting to note that at a particular order there are diagrams which are similar by 
a mere permutaiton of the interaction lines, i.e., by a different time ordering. For example, 
Fig |]b shows four diagrams related by permutations, in the third order of the perturbation 
series. All of these have the same value, and hence the factor of 4 in the rjj term in Eq. |5.2| . 
These permutation factors collaborate with powers of ro in such a way that just the 

same as those for the (Z 2 ) case. There are also diagrams in the third and higher orders (a 
few shown in Fig §d) which correspond to subleading divergences, the removal of which will 
be automatic by their corresponding higher orders. 

We, therefore, see that the j3 function for r has exactly the identical form as that in Eq. 
4.14| for {Z 2 ) c and all the features follow identically. This shows that the phase transition for 
{Z 3 ) c has the same nature as for the (Z 2 ) case. To be more explicit, there exists a transition 
temperature for d > 1 which separates the weak disorder and strong disorder phase for every 
moment. In the high temperature phase (Z 3 ) ~ (Z) 3 and for T < T c , i.e., in the fluctuation 
dominated phase, (Z 3 ) differs from (Z) 3 . This transition temperature is the same for (Z 3 ) 
and (Z 2 ). 

It is now a trivial exercise to extend this for higher moments. |BTJ The effective hamilto- 



nian apart from the free part for 2n chains and mutual 5 function interaction involves the 
following attractive interaction 



r J2 dz <K r 2i-i 2i{z)) 8{r 2 j-i 2j(z)) 

i<3 



Since no new interaction is generated, the f3 function remains the same. [25,30 



VI. RANDOM MULTICRITICAL CASE 



In the previous sections, attention was focussed on the two body interaction case. It is 
known that DPs with pure m-body interaction can also be completely solved. [^,^,0 We 
now investigate the random version of this multicritical case as given by the hamiltonian of 
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Eq. pT3[ As before, we want to evaluate (Z m ) and (ZD- The procedure follows the footsteps 
of the two chain problem, and, therefore, details are skipped. 



A. (Z m ) 



To compute (Z m ), we can perform an averaging over b(z) to obtain, as in section |TJ, an 
m-chain hamiltonian with a pure m-body interaction. The grand universality known for the 
pure system indicates that the multicritical exponents for the binding - unbinding transition 



will be similar to those of Ref [0. For example, for d > 2/ (m — 1), the lengthscale exponent 



would be 2/| e' m |, where e' m = 2 — (m — l)d. ' 

b. (zl) 

A little calculation involving the completion of the square would convince the reader that 
the effective hamiltonian needed for the second moment would involve 2m chains in sets of 
m. It is given by 



n m , m 2 



Wo/ s 2 N , 

1 V 2m d Y .( z )\ m-1 



f dz J2(^-) +v m f dz J] 5[v pp+1 (z)} + 

o i=i\ az J i P =i 



N 2m-l ^ m-1 2m- 1 



V, 



lb _L n f/t_L ^((fcj. 

dz n 5 [ r « g+i ( z )i dz n ^[ r P?+iW] n ^ifH+iwi (s- 1 ) 







5=m+l Q p=l q=m+l 

where f m = v^A and the two sets are represented by p's and g's. The special feature is the 
last term that involves the peculiar m chain - m chain interaction. This generalizes H 2 of 
the two chain case of Eq. [4. 1 b| and |4.1ej . The effect of disorder, so far as the fluctuations 



are concerned, is to introduce a correlation that if m chains meet at a z, the replica would 
also like to enjoy a meeting at that same z. 

The upper critical dimension of r m follows from the dimensional analysis as d m = l/(m — 
1), which is half of the upper critical dimension for the pure case [2/(m — 1) for v m ]. We are 
not sure whether this systematic reduction by a factor of 2 has any deeper significance. 
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For simplicity we choose v m = 0. The perturbation expansion in f m would involve the 
same sets of diagrams as in Fig. 0a except that the propagator now for the thick lines is 
G m (r | z) = (2irz)~^ m ~ 1 ' d ^ 2 m~ d ^ 2 G(r | z/m). With this propagator, the full series can be 
computed. It is transparent to see the occurence of divergences at d — d m . The whole RG 
procedure of section [TV] can be copied in toto by replacing e by e m = 1 — d{m — 1). Hence 
in the multicritical situation we also expect to see a disorder induced phase transition. The 
length scale exponent v m = (2 | e m with an essential singularity for d = d m as in 
Eq. gig . 

VII. SYSTEM WITH MORE THAN TWO CHAINS 

Annealed averaging for the system with two chains with random interaction is simpler 
and not sufficient to give enough information about the effects of disorder. On the otherhand, 
if the above case can be extended to four chains having twobody interaction among each 
other, even the annealed case turns out to be extremely nontrivial. The hamiltonian for the 
four chain system, 

2 

1 r N 4 fdr(z)\ r N 4 
H = 2/ dz Y,{-jt^) +1 dzv (l + b(z)) £ 6(tv(z)) (7.1) 

hj = 1 

i ^ j 

where = v^z) — ?j(z), after averaging, using the Gaussian distribution of b(z), gives the 
following effective hamiltonian 

H cS = y N dzEU {^f) + (¥) 2 + vofo N dzJ2t >j = 1 6(v ij (z)) - 

2v 2 Q AJ N J2 ijjik S{r ij )S(r jk ) - 2v 2 A J N dz^^i 5(r y )«J(r «). (7.2) 

The remarkable feature of the effective hamiltonian is that there are two different kinds of 
attractive interaction one of which involves three chains with a multicritical type interaction 
(Eq. |2.3|) while the other one coupling four chains together, as in the quenched problem. 
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If we take a three chain system, the corresponding effective hamiltonian will involve only 
the three chain term but no four chain interaction of Eq. |7.2| . This term was absent in 
the original three chain hamiltonian. There is now the possibility of a disorder induced 
multicritical behavior, though of pure type. [0-0 

The four chain attractive interaction is marginal at d = 1 and so is the three chain 
interaction. The presence of these two marginal operators is, in general, expected to com- 
plicate the renormalization procedure through their interdependence - but here that does 
not happen. 

The perturbation expansion with the three body and the four body interactions leads to 
three different kinds of diagrams. See Fig [5| The series corresponding to the pure three 



body interaction is already solved. p6| , p!8| The series in the Laplace space involving four 
body interactions (see Fig 5a-c), which contributes to the leading divergence is identical to 
the series for (Z 2 ) c in Eq. |4.4] (Fig. |2]a). The diagrams with mixed three body and four 
body interactions are shown in Fig |5]d,e. In the final series, upto the order shown in Fig |5|, 
the four body and three body contributions get separated into two factors. This shows that 
the resulting renormalization of the two couplings are independent of each other. Because 
of the four body interaction, we expect a disorder induced criticality as for the two chain 
quenched case, but here this happens for a real four chain system - no replica is involved. 
The details and the phase diagram will be published elsewhere. 



VIII. SUMMARY AND DISCUSSION 

We have proposed a random interaction model for two directed polymers and studied the 
first three cumulants of the partition function. We have shown that in the annealed case, 
described by (Z), there can be a disorder induced binding unbinding transition, very similar 
to a pure problem. The exponents are also identical to the pure case. We also pointed out 
certain peculiarities of the annealed problem involving three or four chains. The quenched 
problem is different as reflected through the marginal relevance of the disorder. For d > 1 
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there exists a critical point that demarcates a disorder dominated phase and a pure type 
phase. In the strong disorder phase, there seems to have an extra correlation in the replica 
space, which is absent in the other phase. The lengthscale exponent for the critical point is 
found to be (2 | e where e = 1 — d. It has an exponential divergence at d — 1. Similar 
results were obtained for (Z 3 ). In the replica approach, one needs (Z n ) with n — > which 
does not require interactions other than those which are present in (Z 2 ) and (Z 3 ). Therefore 
the upper critical dimension will remain the same, namely d = 1. 

There are still many open problems, as for example, a replica analysis for this system. 
This requires an explicit expression for (Z n ), correct at least for small n. Such an analysis 
would provide vital information regarding the disorder induced critical point, including 
possible replica symmetry breaking. Is there any other length scale exponent for this new 
crtical point, apart from the one we have calculated? What about other exponents? What, 
if any, is the upper critical dimension of this critical point? A thorough numerical study of 
this system will surely provide valuable insights. 

APPENDIX A: GENERATION OF ATTRACTION: 
PERTURBATION ANALYSIS FOR (Z) AND (Z 2 ) 

In this appendix, we show how the attractive terms in the effective hamiltonian for (Z) 
and (Z 2 ) can be generated perturbatively. 

We proceed to the evaluation of the average of the partition function (Z) by a pertur- 
bation expansion using Eq. |2.1| with the replacement of the short range potential by a 5 
function potential. Formally this leads to the expression for 

Z = J D ri Dr 2 exp (-#„) (1 - H t + H 2 /2\ + ...) (Al) 

where Hq corresponds to the free part of the two chains, and Hi represents the interaction 
part, J dz vq (1 + b(z)) 5 (r 12(2)). A little manipulation after substituting the explicit form 
of Hi gives the connected partition function as 
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Z c = J* dz x v (1 + b(zt)) J {ry} G(r N - n | N - z x ) G(r 1 - r | z x ) x 
G(r' N -rtlN-zt) G(n - r' | z x ) + 
J N dzt / Z1 dz 2 vl (1 + b{ Zl )){l + 6(22)) x 
I{ r ,r'} G ( r N ~ ri I N - z x ) G(r 12 I ^12) G(r 2 - r | z 2 ) x 
G(rV - n I iV - zi)G(na | z 12 )G(r 2 - r' | 22) + - (A2) 

where G(r | z) = (2nz)^ d ^ exp(— r 2 /2z) is the distribution function for a (^-dimensional 
Gaussian chain of chain length z and end-end distance vector r and J^ rr ij the integrations 
over all dummy spatial coordinates. The convention is to use r and r' for the two chains 
and Zij = Zi — Zj. The factorials in the denominators of the terms of Eq. [XT] are absorbed 
by introducing "time" (z) ordering, which restricts Zi + \ < Zi, in the integrals. Diagrams 
upto second order corresponding to the series in Eq. [A^ are shown in Fig ^. The familiar 
normalization / drG(r | z) = 1 and the integrations over the spatial end coordinates lead to 
much more simplification for Z c . The explicit form of Zc upto two loop term is given by 

Z c = V J N dz lVo (l + b{z x )) - J N dz x ; 21 dz 2 J {r} vl (1 + b( Zl ))(l + 6(22)) x 
G 2 (r 12 I z n ) + dzi Jo* dz 2 J Z2 dz 3 J {r} vl(l + b( Zl ))(l + b(z 2 )) x 

(1 + 6(z 3 ))G 2 (r 12 I z 12 )G 2 (v 23 \ z 23 ).... (A3) 

where V is the transverse volume. The actual meaningful quantity (Z c ) is computed from 
Eq. |A3| after averaging it with the distribution P(b) of Eq. |2.2a| . This gaussian distribution 



with zero mean ensures that any term involving odd number of 6(z)'s should vanish after 
averaging. Therefore, the contribution from the first order term, VqNV, is only from the 
pure part. In the one loop level of Eq. IA3], there are two surviving terms after disorder 



averaging. One is the pure term which does not require any averaging and its contribution 
to (Z c ) is 

- v2 ° v{An) ~ d/2 Y^W) N2 ~ d/2 - (A4) 

The other nonvanishing part which contains an even number of disorder interaction is 
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— Uq / dzi ( dz 2 [ dr 1 dr 2 b(zi) b(z 2 ) G 2 (r 12 | 212). 
Jo Jo J 

After averaging, the two points Z\ and z 2 along the chain merge together giving rise to a 
term 

rN 

— / dzi vl A5 2 (ri - r 2 ) dr 1 dr 2 (A5) 
Jo 

where we used the fact that G(r | 0) = 5(r). Because of this merging of the two points 
along the chain, this term contributes to the first order term but with a negative sign which 
shows the presence of a newly generated attraction. In other words, a second order term for 
a particular realization (before averaging), looks like an attractive first order term after DA. 



Proceeding in the same fashion we can evaluate the two loop term of Eq. |A3| . In the 
two loop part, the nonvanishing contributions are: 

rN rZ\ I-Z2 



W ^0 



fjy r z i r z 2 r 

/ dzi / dz 2 / dz 3 / G 2 (r 12 | z 12 )G 2 (r 23 \ z 23 ). 

JO Jo JO J {r} 



This term involving only pure type interaction, after integration over spatial coordinates, 
gives 

rN rz\ rz 2 



ri\ rz\ rZ2 

V^(47r)- d / dz x / dz 2 / dz, 
Jo Jo Jo 



^ Zl - Z2 y d i\z 2 - Z3 )- d / 2 



(A6) 



fN rzi rZ2 r 

(ii)vl / dzt / dz 2 / dz 3 / b(z 1 )b(z 2 ) G 2 (r 12 \ z 12 )G 2 (r 23 \ z 23 ). 

Jo Jo Jo J {r} 

As in the earlier case this reduces to a one loop term after averaging as 

v 3 A (An)- d / 2 [ N dz x f 1 dz 2 r dz 3 I 5 2 {v 12 ) (z 2 - z 3 )- d ' 2 . (A7) 

Jo JO JO J{r} 

(iii) The third nonvanishing contribution is from 

rN rZi rZ2 r 

v 3 / dzi / dz 2 / dz 3 / b{z 2 ) b{z 3 ) G 2 (r 12 I z 12 ) G 2 (r 23 \ z 23 ). (A8) 

Jo JO JO J{r} 

This term after averaging becomes 
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v 3 A (47r)- rf / 2 r d Zl f 1 dz 2 [ 5 2 (r 23 ) { Zl - z 2 )- d ' 2 . (A9) 

JO JO J{r} 

There is one more term involving b(zi)b(z 3 ) which vanishes after averaging because of the 
specific time ordering which rules out the merging of z\ and 23. It can be easily be checked 
that this merging of two "random interaction" lines (wavy lines) into one single pure line 
and the subsequent reduction of order occurs at each order (> 1) involving consecutive 
pairs of even number of "random" lines. Thus an attraction is generated at each order very 
systematically. The coupling constant of this term is proportional to t>o 2 A. In some of the 
above expressions, [ |A5|, |A7|, [A9|| the presence of the 5 2 term needs special attention since 



it is illdefined even in the theory of distribution sense. We can avoid this problem by 
taking a spread out 5 function, and then taking the limit at the end. This would change 
the coupling to VqA/Q, where Q is the arbitrary "spread out" or cutoff volume. Since one 
gets back a single 6 function, it can be associated with the pure term, thereby changing the 
problem to a pure one with a reduced coupling constant v = v — VqA/Q. See Eq. P^2] . 
Another way to tackle this difficulty is to start with a short range potential V(r), and appeal 
to RG arguments as done in section fH\ . 

From now onwards it's time for (Z 2 ) A few diagrams for (Z 2 ) with vq = are shown in 
fig. 7a. In the first order, the only diagram which has nonzero contribution to (Z 2 ) c is Fig. 
7al. This contribution after disorder averaging is V 2 VqAN. As was mentioned in the text 
this V 2 has come from the independent spatial integrals. Proceeding in the similar fashion, 
we write 

figlai = v / dz\ [ dz 2 b(zi)b(z 2 ) [ G 2 (r 12 | z 12 ) x 

Jo Jo J{r} 

/ dz[ dz' 2 b(z[)b(z' 2 ) G 2 (r' 12 \z' 12 ). 

JO Jo J{r,r'} 

In disorder averaging, the only relevant contribution comes from the pairing of b(zi) b(z[) 
and b(z 2 ) b(z 2 )- The other possibility in which b{z\) b(z 2 ) and b(z[) b(z' 2 ) are paired up (Fig 
7b2) is not considered here since this generates Vq type terms which are not to be included 
for the ^0 = case. After appropriate disorder averaging the above expression becomes 
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^A 2 V 2 4vr d r(e)s- (2+£) , (A10) 

where e = 1 — d. 

To make the evaluation, after DA, easier we follow a different convention for the diagrams, 
Fig 0b and 2a. The thick line represents the two members of a pair jointly and is represented 
by G 2 (ri — r 2 \ z\ — z 2 ) for the ends ri and r 2 at which the two chains are tied at lengths 
z\ and z 2 . For example, the diagram which corresponds to the last expression (Eq. |A10| ) is 
given by Fig [7]b3, which is also Fig 2a2. (Note fo = v 2 A.) 

The next diagram which is important in the next higher order is given in Fig 7a4, the 
contribution of which can equivalently be calculated from Fig 7b5 or 2a3 as 

rN pZl pZ'2 p 

(v 2 A) 3 / dzi / dz 2 / dz 3 / G 2 (r 12 | z 12 )G 2 (r 23 \ z 23 ) x 

JO JO Jo J{r, r'} 

G 2 (r' 12 I z 12 )G 2 (r' 23 \ z 23 ). (All) 

In the Laplace space this becomes 

K 2 A) 3 V 2 r 2 (e)(47rr 2 V( 2+2e ) (A12) 

The diagrams having odd numbers of wiggly lines trivially vanishes after DA. This can be 
generalized to arbitrary orders since only ladder type diagrams are involved. Eq. [14] would 
follow by substitution f = v 2 A. 

APPENDIX B: DRESSED PROPAGATORS AND {Z 2 ) 

We first show the two different dressed propagators. The one for which both the chains, 

tied at the ends (r,z) and (r',z'), meet each other n times at (ri,zi), {r 2 ,z 2 ) 1 (r n , z n ) 

(Fig |]d) is given by 

G$(r -v'\z-z') = v% J {r} f z z , d Zl \ Z J dz 2 ... J;,"" 1 dz n G 2 (r - Tl \z-z 1 )x 
G 2 (r 12 I 2i 2 )....G 2 (r n _i n | z n -i n )G 2 (r 

n r z n Z I 

Use of the identity G 2 (r | z) = (4ivz)~ d / 2 G(r \ z/2) and the Markovian property 
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J rfr 2 G(ri - r 2 | z 1 )G{r 2 - r 3 | z 2 ) = G(ri - r 3 | z x + z 2 ) 

leads to the following expression for G^(r — r'\z — z'). 

G { M ] (r -r'\z-z') = v^- nd ' 2 G{r - r' \ (z - z')/2) x 
J* d Zl & dz 2 ... r z r' dz n (z - Zl )- d / 2 z u d/2 .... x 

-d/2 , _ i\-d/2 

So we need the following integral 

f d Zl (z - Zl )- d ' 2 p dz 2 ( Zl - z 2 )- d/2 ... dz n {z n . x - z n y d ' 2 (z n - z')- d ' 2 . 

Jz' Jz' J z' 

A change of variable 

Z{ Z{ z 

and use of the convolution theorem in the Laplace space straightaway yields 

s - e ' ( „+i) r „+i^ 

where s is the Laplace conjugate to the chainlength (more precisely to z = z—z'). Converting 
this to inverse Laplace space and combining all other factors, the final form of such a 
propagator becomes 

Gi\r | z) = (-v r(4n)-^ d / 2 z^y-i / G(r \ z/2). (Bl) 

T[{n+ lje'j 

Following the same track, the other dressed propagator, for which the two member chains 
are tied only at one end say (r, z) other than n meetings, has the form (see Fig 2c) 

g ( o\z) = (-^n^r^'H^ ^+q ( B 2) 

This propagator is independent of any space coordinate because of the spatial integration 
over the open end coordinates. 

To take care of arbitrary number of meetings we sum over n. Hence the final dressed 
propagators are 
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oo 



n=0 



oo 



G M (r - r' \ z - z') = J2 G 

n=0 

r2\ 



in) i i | / 

)y (r — r | z — z 



Therefore the series for (Z ) c can be written as 

/•JV 

(Z 2 ) c = f / d^G^iV - Zl ) G 2 ( Zl )+ (B3) 
fjj / d2j f* cfe 2 / G%(iV - 2 X ) G M (r 12 | z 12 ) G M (r' 12 | 2 12 ) G%{z 2 ) + .... (B4) 

JO JO J {r, r'} 



In fact, it is possible to write the whole series to all orders in perturbation. Completing the 
integrations over the spatial coordinates we get 

(Z 2 ) c = V 2 f J" d Zl G 2 (N - z 1 )G 2 (z 1 ) + f 2 V 2 f N d Zl ; 21 dz 2 G 2 {N - Zl ) x 
Gl(z 12 )Gl(z 2 ) + f 3 V 2 f N dz x Jo 1 dz 2 ^ dz 3 Gl(N - Zl ) x 

G 2 m (z 12 )G 2 m (z 23 )Gq(z 3 ). (B5) 

Here Gm(z) is obtained from G«(r | z) after integration over the spatial coordinate. We 
use the convolution theorem for Laplace transforms which leads to the following expression 
for Z the Laplace transform of (Z 2 ) c (Eq. f4.2| ) 

Z |^o= V 2 [f £ 2 ( S ) +f 2 £ (s)£Af(s)£o(s) + 

flQ {s)Q 2 M {s)Q {s) + ....]. (B6) 

where, Z = / °° e - sAr (Z 2 ) c , £ p = / °° e~ sN G 2 p (z), with p being O or M. The Laplace trans- 
forms Qo(s) and are given by 

v (4 7 r)-^+^) rf / 2 < i+ " 2 r^+^(6Qr((n 1 + n 2 y + 1) 

Wj n f^ 2 r(l + me')r(l + n 2 e') s (m+n 2 )(i-d/2)+i ^ 

and 

(4 7r) -(n 1+ n 2 +2)d/2^n 1 +n 2r n 1+ n 2(e / )r[(Tli + + £ ] 

[S) nb> + ni)e']r[(l + n 2 )e'] sC^y+e " 

Substituting these expressions in equation |B^, we get back equation for Z \v ^o in the 
Laplace space. The results for ^0 = can be obtained from the first term of each sum, i.e., 
for Tlx = n 2 = 0. 
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FIGURES 

FIG. 1. Flow diagrams for coupling constant u in different dimensions. u*(= Aire') represents 
the nontrivial fixed point. 

FIG. 2. (a) The only contributing diagrams in {Z 2 ) \v =o upto third order. Only ladder dia- 
grams occur. A thick line corresponds to a pair of chains. A wiggly line stands for an fo factor in 
the evaluation of the diagrams, (b) A typical diagram for (Z 2 ) \v ^o- The dots on the thick lines 
represent intrapair interactions (fo)- (c) The dressed propagator with two chains tied at only one 
end. The dashed lines represent the mutual (5-function type interaction with coupling constant vq. 
(d) The dressed propagator with two chains tied at both (r, z), and (r', z'). 

FIG. 3. Flow diagram for r in various dimensions. r*(= — e) represents the nontrivial unstable 
fixed point. For the mth order multicritical case r is to be replaced by r m . The three figures would 
be for d < d m , d = d m , and d > d m , where d rn = l/(m — 1). 

FIG. 4. The second (a), third (b) and fourth (c) order (in fo) connected diagrams for (Z 3 ) \v =o- 
In fourth order, there are a few other similar diagrams which contribute to the leading divergence. 
For connectedness, the series has to start at order two. (d) A few diagrams which contribute to 
subleading divergences in the third and fourth order in fo 

FIG. 5. Four chain diagrams for the annealed problem with four chains. The wiggly lines 
represent fo type interaction, and a solid horizontal line connecting three chains is the three chain 
^-function interaction. The first three terms of the series involving only four chain interactions are 
shown in (a), (b), and (c). (d and e) Two cases involving both the three and four chain interactions, 
(f) A possible diagram in second order with different chain combinations for the interactions. This 
contributes in the subleading divergence. 

FIG. 6. Diagrams for (Z) upto second order in fo and vob(z). The wavy and dotted lines 
interactions with coupling constants vob(z) and vq respectively. 



32 



FIG. 7. (a) Diagrams involving only vob(z) for {Z 2 ). (b) Diagrams of (a) after disorder aver- 
aging. (See the caption of Fig ||.) Fig (a2) with odd number of wavy lines vanishes after DA. 
Different pairings lead to two possiblities (b2) and (b3) from (a3). Similarly for a4, there are two 
diagrams (b4) and (b5) after DA. Diagrams (b2) and (b4) are not considered for the (Z 2 ) \v =o 
case. 
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